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Abstract 

The determination of renormalization factors is of crucial importance in lattice 
QCD. They relate the observables obtained on the lattice to their measured counter- 
parts in the continuum in a suitable renormalization scheme. Therefore, they have 
to be computed as precisely as possible. A widely used approach is the nonper- 
turbative Rome-Southampton method. It requires, however, a careful treatment of 
lattice artifacts. In this paper we investigate a method to suppress these artifacts 
by subtracting one-loop contributions to renormalization factors calculated in lattice 
perturbation theory. We compare results obtained from a complete one-loop sub- 
traction with those calculated for a subtraction of contributions proportional to the 
square of the lattice spacing. 

1 Introduction 

Renormalization factors in lattice Quantum Chromodynamics (QCD) relate observables 
computed on finite lattices to their continuum counterparts in specific renormalization 
schemes. Therefore, their determination should be as precise as possible in order to al- 
low for a reliable comparison with experimental results. One approach is based on lattice 
perturbation theory [I]. However, it suffers from its intrinsic complexity, slow conver- 
gence and the impossibility to handle mixing with lower-dimensional operators. Therefore, 
nonperturbative methods have been developed and applied. Among them the so-called 
Rome-Southampton method [2] (utilizing the RI-MOM scheme) is widely used because of 
its simple implementation. It requires, however, gauge fixing. 

Like (almost) all quantities evaluated in lattice QCD also renormalization factors suffer 
from discretization effects. One can attempt to cope with these lattice artifacts by extrap- 
olating the nonperturbative scale dependence to the continuum (see Ref. [3]) or one can 



try to suppress them by a subtraction procedure based on perturbation theory. Here we 
shall deal with the latter approach. 

In a recent paper of the QCDSF/UKQCD collaboration [I] a comprehensive discussion 
and comparison of perturbative and nonperturbative renormalization have been given. 
Particular emphasis was placed on the perturbative subtraction of the unavoidable lattice 
artifacts. For simple operators this can be done in one-loop order completely by comput- 
ing the corresponding diagrams for finite lattice spacing numerically. While being very 
effective this procedure is rather involved and not suited as a general method for more 
complex operators, especially for operators with more than one covariant derivative, and 
complicated lattice actions. An alternative approach can be based on the subtraction of 
one-loop terms of order a 2 with a being the lattice spacing. The computation of those 
terms has been developed by the authors of Ref. [6j and applied to various operators for 
different actions. In this paper we use some of those results for the analysis of Monte Carlo 
data for renormalization coefficients. 

We study the flavor-nonsinglet quark-antiquark operators given in Table [TJ The cor- 
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Table 1: Operators and their representations as investigated in the present paper. The 
symbol {...} means total symmetrization. A detailed group theoretical discussion is given 
in 0. 

responding renormalization factors have been measured (and chirally extrapolated) at 
(3 = 5.20, 5.25, 5.29 and 5.40 using Nf = 2 clover improved Wilson fermions with pla- 
quette gauge action [I]. All results are computed in Landau gauge. The clover parameter 
csw used in the perturbative calculation discussed below is set to its lowest order value 
c sw = 1- 

2 Renormalization group invariant operators 

We define the renormalization constant Z of an operator O from its amputated Green 
function (or vertex function) T(p), where p is the external momentum and the operator 
is taken at vanishing momentum. The corresponding renormalized vertex function and 
the Born term (with all lattice artifacts included) are denoted by T R (p) and r Born (p), 
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respectively. If there is no mixing, Z can then be obtained by imposing the condition 



Itr^r^tp)- 1 ] = 1 (1) 

for vanishing quark mass at p 2 = /i 2 , where fi is the renormalization scale. The Z factor 
relates the renormalized and the unrenormalized vertex function through 

T R (p) = Z^ZT(p), (2) 

with Z q being the quark field renormalization constant determined by 

7 / a _ tr R Ea 7a sin(ap A ) aS" 1 ^)] 

in the chiral limit again at p 2 = /i 2 . Condition ([1]) together with ([3]) defines the RI'-MOM 
renormalization scheme. Here S 1-1 is the inverse quark propagator. Using ([TJ) we compute 
Z from 

Z; 1 Z±tr[T(p)T B °™(p)- 1 ]=l. (4) 

For operators transforming as singlets under the hypercubic group H(A), such as O s , Z 
can depend on the components of p only through H (4) invariants. 

For operators belonging to an H (4) multiplet of dimension greater than 1 the condition 
([TJ) violates H(A) covariance and would in general lead to different Z factors for each 
member of the multiplet. In Ref. |1] an averaging procedure has been proposed to calculate 
one common Z factor for every multiplet. Labeling the chosen operator basis by i = 
1,2, ... ,d the common Z was calculated from 

z ? z \T,^ tT N^rfnp)- 1 ] = i • (5) 

i=l 

This condition leads to an if (4)-invariant Z for the operators without derivatives in Table [TJ 
However, in general this is not the case. 

It is not difficult to devise a renormalization condition that respects the hypercubic 
symmetry. Choosing a basis of operators (again labeled by i), transforming according to a 
unitary irreducible representation of H(4), the relation 

E?,itr[rf "(p)rf™(p)t] 

defines a Z factor which is invariant under i?(4), provided that the quark field renormal- 
ization factor is also H(4) invariant. The derivation of renormalization condition ([6]) is 
given in the Appendix. For the operators without derivatives the definitions fl6]) and (j5]) 
are equivalent. For the considered operators with one derivative the resulting differences 
turn out to be negligible. In the following the Z factors will be determined from ()6]) using 
the operator bases given in Table [TJ This is our version of the RI'-MOM scheme. 
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We define a so-called RGI (renormalization group invariant) operator, which is inde- 
pendent of scale M and scheme S, by [I] 

C RGI = AZ S (M) O s (M) = Z KG \a) £> bare (7) 

with 

and the RGI renormalization constant (depending on a via the lattice coupling) 

Z* G \a) = AZ s (M)Zt re (M,a). (9) 

Here g s , •y s and j3 s are the coupling constant, the anomalous dimension and the (3- function 
in scheme S, respectively. Relations (J7J), (jSJ) and (jUJ) allow us to compute the operator O 
in any scheme and at any scale we like, once Z RGl is known. Therefore, the knowledge 
of Z RG1 is very useful for the renormalization procedure in general. Ideally, Z RGl depends 
only on the bare lattice coupling, but not on the momentum p. Computed on a lattice, 
however, it suffers from lattice artifacts, e.g., it contains contributions proportional to a 2 p 2 , 
(a 2 p 2 ) 2 etc. For a precise determination it is essential to have these discretization errors 
under control. 

As the RI'-MOM scheme is in general not 0(4)-covariant even in the continuum limit, 
it is not very suitable for computing the anomalous dimensions needed in (JBJ). Therefore 
we use an intermediate scheme S with known anomalous dimensions and calculate Z 1 
as follows: 

Z RG \a) = AZ S (M = /,) Z^ MOM (M = /,) z£.- MQM (/i, a) . (10) 

It turns out that a type of momentum subtraction scheme is a good choice for S (for 
details see Ref. 0]). The formula which is used to compute the transformation factor 
^ri'-mom (/•*) i s gi ven there together with all needed coefficients of the /3-function and 
anomalous dimensions, which are based on continuum three-loop calculations such as those 

in [HE Eg. 

On a lattice with linear extent L the scale fi should ideally fulfill the relation 

l/L 2 «A^ CD «/! 2 «l/a 2 . (11) 

In that case Z RG1 (a) would be independent of /i, and from the resulting plateau we could 
read off the corresponding final value. However, in practice a\i is not necessarily small 
leading to non-negligible lattice artifacts that have to be tamed. A promising tool to 
control lattice artifacts in a systematic way is lattice perturbation theory: We expect that 
after subtracting these perturbative terms the calculation of the Z factors can be done 
more accurately. 
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3 Subtraction of all lattice artifacts in one-loop order 



In standard lattice perturbation theory the one-loop renormalization constants are given 
in the form 

Z(^a) = l + ^-% (7omM+A), C F = \ . (12) 

16 7T Z 3 

This means that the a-dependence is retained only in the logarithm and implicitly in g, 
while in all other contributions the limit a — > has been taken. 

However, there is no need to do so. We can keep a finite everywhere and thus evaluate 
the lattice artifacts at one-loop order completely, proceeding as follows. Let us denote by 
F(p, a) the total one-loop correction and by F(p, a) the expression resulting from F(p, a) 
by neglecting all contributions which vanish for a — > 0. The difference 

D(p,a)=F(p,a)-F(p,a) (13) 

represents the lattice artifacts in one-loop perturbation theory and is used to correct for 
the discretization errors: 

ZZ: MOM (P, «)MC, sub = ZZ; M ° M (P, a) MC - ^- 2 C F D(p, a) . (14) 

There is a certain freedom in choosing the coupling in (j!4j) . It turned out that the use 
of the boosted coupling 

»S = = 9 2 + 0{g A ) (15) 

(P(g) being the measured plaquette at (3 = 6/g 2 ) is quite successful in estimating the 
higher-order discretization effects. With the prescription (|14p all lattice artifacts in one- 
loop order are subtracted. 

In Fig. [1] we show the effect of subtraction on the RGI renormalization factors for 
selected operators of Table [H For all operators we recognize after subtraction a remarkable 
smoothing and a pronounced plateau as a function of p 2 for p 2 > lOGeV 2 . The large 
bending in the small p 2 region might indicate the breakdown of perturbation theory (cf. the 
discussion in [3]). The examples show that the one-loop subtraction of lattice artifacts (I14p 
works very well and, moreover, is needed for a precise determination of the renormalization 
constants. The final values for Z RGI from (111) I) are obtained by a fit with an ansatz [4J 

2&-mom(p) Z^ MOM (p, a) M c, sub = Xz^m^kWf] + Cl aV • (16) 

The free parameter b\ takes into account that the transformation factor -^ri/_mom(p) ^ s 
known to three- loop order (g 5 ) 6 only. Further possible lattice artifacts are parametrized 
by ci a 2 p 2 . 

For practical reasons the numerical calculation of F(p, a) - and therefore the calculation 
of Z RGl using f lT6|) - is restricted to operators with at most one derivative and for Nf = 2 
only. In order to perform the subtraction for a wider class of operators and/or for Nf = 2+1 
(where the considered lattice action becomes more complicated) we have to look for an 
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Figure 1: Z^ G1 (left) and Z^ Gl (right) for (3 = 5.40. The Z factors obtained without 
subtraction are shown as red squares, those with complete one-loop subtraction (fT4"l) as 
blue triangles. (The necessary scale transformation factors for the momenta are given at 
the end of Section HI) 



alternative method. One possibility which will be discussed in the next sections is a 
'reduced' subtraction: Instead of subtracting the complete one-loop lattice artifacts we 
subtract only the one-loop terms proportional to a 2 , if they are known for the given action. 



4 Subtraction of order a 2 one-loop lattice artifacts 
4.1 Lattice perturbation theory up to order g 2 a 2 

The diagrammatic approach to compute the one-loop a 2 terms for the Z factors of local 
and one-link operators has been developed by some of us [HI ED]- The general case of 
Wilson type improved fermions is discussed in [IT] . For details of the computations we 
refer to these references. Here we give explicitly the results for the operators and actions 
investigated in this paper (massless improved Wilson fermions with c$w = 1, plaquette 
gauge action, Landau gauge). 

Using the relation (|6]) we compute a common Z factor for each multiplet given in 
Table [TJ The results £1X6 £IS follows: 
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Here we have introduced the notation 



Px 



A=l 



with p\ being the momentum components. Note that terms of type (S^/ S 2 )\og(a 2 S 2 ), 
appearing in Zg, Zy, Za, Zt, all have the same coefficient which arises solely from the 
quark wave function renormalization constant Z q . The corresponding one-loop vertex 
functions Tj(p) in (jBJ) do not contain such a structure. For later purposes we write the 
Z factors generically as 



1 . 9 2 C-F ry 2 2 rv{a 2 ) i s 



(19) 



We emphasize that the numerical coefficients in the above expressions are either exact 
rationals or can be computed to a very high precision. 

In Figs. El [3] and [4] we present a 2 g 2 Z l f_ i oop (p, a) for selected operators as a function of 
a 2 p 2 on a finite lattice, where we choose the lattice momenta as p\ = (2iti\)/(aL\). Here, 
i\ are integers and L\ is the lattice extension in direction A. We compare the correction 
terms for a general set of momenta with those obtained for the momenta used in this 
investigation at f3 = 5.40 on 24 3 x 48 lattices and with 'diagonal' momenta, i.e., momenta 
on the diagonal of the Brillouin zone. 



7 




Figure 2: a 2 g 2 Z[°i^ oop (p, a) for operators O s (left) and O v (right) as a function of a 2 p 2 on 
a 24 3 x 48 lattice at /3 = 5.40. The green filled circles are the values for an arbitrary set 
of (mostly non-diagonal) momenta, whereas the red filled squares are obtained from the 
momenta used in this investigation. The blue line is computed from diagonal momenta. 
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The figures show that the momenta of the actually measured Z factors are very close to 
the diagonal. Furthermore, one recognizes that the magnitude of the calculated one-loop 
a 2 corrections in the used momentum range is small but not negligible compared to the 
measured values which are of order 1 (see also Fig. [T]). Therefore, one can expect that the 
subtraction of those terms yields a noticeable effect. 



4.2 Subtraction of lattice artifacts up to order a 2 

The subtraction procedure of order a 2 terms is not unique - we can use different definitions. 
The only restriction is that at one-loop order they should agree (treating Z^ I ~ MOM (p, o)mc 
in perturbation theory). We investigate the following possibilities, 
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where g+ can be chosen to be either the bare lattice coupling g or the boosted coupling 
#b ( TT5T) . (In the following we denote subtraction type ( |20l) by (s) and ( |2TI) by (m)). 
With ansatz (s) the one-loop a 2 correction is subtracted 'directly' from Z^l~ M M (p, o)mc- 
Subtraction type (m) factorizes the one-loop a 2 correction from the nonperturbative Z 
factor. 

The Z RGl are computed from (fTUj) using (s) or (m), where we expect slightly dif- 
ferent numbers depending on the choice of coupling g*. The only significant errors to 



?RI'-MOM 
^bare 



(p, a)Mc,sub are due to the Monte Carlo simulations. 



In Fig. we show how the subtraction of lattice artifacts (complete and a 2 ) affects 
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Figure 5: Unsubtracted and subtracted renormalization constants for the scalar operator 
O s (left) and the tensor operator O t (right) at f3 — 5.40, for p 2 > 10 GeV 2 and r Aj^g = 
0.700. The complete subtraction is based on ( 1T4"|) . whereas the a 2 subtractions are of type 
(s) and (m) with g+ = g B . 



the renormalization constants for the scalar and tensor operators. The complete one-loop 
subtraction results in a clear plateau for both Z RG1 factors. Using the a 2 subtractions 
there remains a more or less pronounced curvature which has to be fitted. From the 
definitions of the subtraction terms it is clear that they vanish at a 2 p 2 = 0. Moreover, for 
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small p 2 ~ 10 GeV 2 the subtraction methods (s) and (JHJ) already agree, as they should. 
However, as discussed above, Z RGI can only be determined from sufficiently large momenta 
(p 2 ^ 10 GeV 2 ), where differences arise between the various procedures. Therefore the 
results for Z 1 may differ depending on the kind of subtraction. As can be seen in Fig. [5] 
this effect varies strongly from operator to operator. 

4.3 Fit procedure 

Compared to the complete one-loop subtraction we expect that Z R l r ~ MOM (p, a)Mc,sub as 
computed from (s) or (m) contains terms proportional to a 2n (n > 2) even at order g 2 , as 
well as the lattice artifacts from higher orders in perturbation theory, constrained only by 
hypercubic symmetry. Therefore, we parametrize the subtracted data for each (3 in terms 
of the hypercubic invariants S n defined in ( ITS]) as follows 

s RI' mom Z RGl (a) 

^RI'-MOmCp) Z bare (P> a )MC,sub = [l_ + ^ gS)8l + ^ 

a 2 (ci S 2 + c 2 |p + c 3 + fl4 ( c 4 (S 2 ) 2 + c 5 S4) + a 6 (c 6 (5 2 ) 3 + c 7 S 2 + c 8 5 6 ) . 

There are also further non-polynomial invariants at order a 4 , a 6 , but their behavior is 
expected to be well described by the invariants which have been included already. Ansatz 
( 122]) is a generalization of ([IB]) : After the 'reduced' one- loop subtraction of lattice artifacts 
the Z factors are expected to depend more strongly on a 4 or a 6 hypercubic invariants than 
after the complete one-loop subtraction (see Fig. [5]). The parameters Ci, . . . ,c 8 describe 
the lattice artifacts. 

Together with the target parameter Z RGl (a) we have ten parameters for this general 
case. In view of the limited number of data points for each single (3 value (5.20, 5.25, 5.29, 
5.40) we apply the ansatz (122]) to several (3 values simultaneously with 

Z RG \a) ^ Z^\a k ) 

AZS(p) [1 + h (^)8] AZ s (p) [i + bl ( ^)8] > ^ ) 

where k labels the corresponding (3 value (a^ = a(/3fc)). The parameters q are taken to 
be independent of (3. This enhances the ratio (number of data points)/(number of fit 
parameters) significantly and we obtain several Z R (%) at once. The fit is performed 
by a nonlinear model fit which uses - depending on the actual convergence - either the 
Nelder-Mead or a differential evolution algorithm [12] . Additionally, we have checked some 
of the fit results using MINUIT p3]. 

The renormalization factors are influenced by the choice for r A^. This quantity 
enters AZ S (M) in §E§ via the corresponding coupling g s (M) (for details see [1]). We choose 
tq A^jg = 0.700 [H]. In order to estimate the influence of the choice of r A^jg we also use 
r Af^g = 0.789 calculated in [15]. The Sommer scale r is chosen to be r = 0.501 fm and 
the relation between the lattice spacing a and the inverse lattice coupling (3 is given by 
r /a = 6.050 ((3 = 5.20), 6.603 ((3 = 5.25), 7.004 ((3 = 5.29) and 8.285 ((3 = 5.40) [16]. 
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5 Renormalization factors for local and one-link op- 
erators 



The fit procedure as sketched above has quite a few degrees of freedom and it is essential to 
investigate their influence carefully. A criterion for the choice of the minimal value of p 2 is 
provided by the breakdown of perturbation theory at small momenta. The data suggest [4J 
that we are on the 'safe side' when choosing p^ in = 10 GeV 2 . As the upper end of the fit 
interval we take the maximal available momentum at given coupling (3. 
Other important factors are 

• Type of subtraction: As discussed above the procedure of the one-loop subtrac- 
tion is not unique. We choose different definitions (s) and (m) with either bare g or 
boosted coupling 

• Selection of hypercubic invariants: For the quality of the fit it is essential how 
well we describe the lattice artifacts which remain after subtraction [TT1 [T8] . This is 
connected to the question whether the a 2 subtraction has been sufficient to subtract 
(almost) all a 2 artifacts. Therefore, we perform fits with various combinations of 
structures with coefficients q in (12"2l) . One should mention that the concrete optimal 
(i.e. minimal) set of q depends strongly on the momenta of the available Monte 
Carlo data - nearly diagonal momenta require fewer structures to be fitted than far 
off-diagonal ones. 

The analysis should provide an optimal restricted set of parameters which can be used 
as a guideline for other classes of operators. Nevertheless, one has to inspect every new 
case carefully. 

The results for Z RGI will depend on the above mentioned factors. As a detailed pre- 
sentation for all operators and /3-values would be too lengthy, we select some operators 
and/or /3 values and take the corresponding results as a kind of reference. All results pre- 
sented in this section are computed for r Aj^g = 0.700. The choice r Ajjg = 0.789 leads 
to qualitatively similar results. The large number of parameters in ansatz (j22|) calls for a 
combined use of the data sets at (3 = (5.20, 5.25, 5.29, 5.40) for our fit analysis as indicated 
in ([25]) . With the choice p^ in = 10 GeV 2 this results in 94 data points available for the 
corresponding fits. Additionally, we should note that the errors on our fit parameters are 
those obtained from the nonlinear model fit. They differ from the error calculation for the 
Z RGI based on ( ITS]) and used in [4]. 

5.1 Dependence on the subtraction type 

In Fig. |6]we present the Z RGI for operators O s , O , O t and V2 ' a for the different subtrac- 
tion types using the fit ansatz ( f2~2l) with all q 7^ 0, i.e., we include a 2 , a 4 and a 6 terms. From 
the discussion in Section 14721 we expect that the resulting differences vary from operator to 
operator (cf. Fig. [5]). 

From Fig. [6] we observe that the complete one-loop subtraction (1) and the subtraction 
(2) agree within 1 %. This is not unexpected because the subtraction schemes are similar 
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Figure 6: Z RG1 of selected operators at /3 — 5.40 as a function of the subtraction type 
(subtype): 1: complete subtraction (114"]) with g± = gs, 2: (s) with = 3: (m) with 
— Psj 4: (s) with = g>, 5: (m) with g± = g. The horizontal borders of the shaded 
area show a 1% deviation from case 1. 



and the gauge couplings coincide. The differences in the results for (2) and (3) can be 
used as an indication for a systematic uncertainty in the determination of Z RG1 based 
on the schemes (s,m). We observe that both subtraction approaches are numerically 
almost equivalent. Choices (4) and (5) lead to Z 1 factors which are partly outside the 
1 % deviation. Generally, we recognize that all subtraction procedures for both bare and 
boosted couplings produce fit results within a reasonable error band width. 

In order to test the effect of subtraction we compare the g 2 a 2 contributions as given in 
( I17p with the remaining lattice artifacts of the Monte Carlo data fitted after subtraction, i.e. 



the result for (J22J) setting Z 



RGI/ 



0. In Fig. [7] we show the results for the same selected 



operators choosing gs- In the small p 2 region the remaining lattice artifacts are significantly 
smaller than the one- loop a 2 terms (operators O s , O t and V2 ' a ). In case of already small 
one-loop a 2 artifacts (operator O v ) the final artifacts remain small. This behavior strongly 
suggests to subtract the one-loop a 2 terms before applying the fit procedure. 

Since the boosted coupling g B is assumed to remove large lattice artifacts due to tadpole 
contributions in the perturbative series, we will use gs in the following. In addition, we 
restrict ourselves to subtraction type (s), which is closest in spirit to the complete one- loop 
subtraction studied in [4] (leading approximately to a plateau in the Z RG1 as a function of 
p 2 ). 
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Figure 7: Lattice artifacts for Z RG1 of selected operators for /3 = 5.40 as a function of p 2 
choosing g* = gs- The blue filled circles are the corresponding g 2 a 2 correction terms, the 
red open circles are the fit results for (|22|) setting Z RG1 (a) = 0. 



5.2 Dependence on hypercubic invariants 

Now we discuss the dependence on the hypercubic invariants included in the fit ansatz 
( 122|) . The goal is to select a reasonable set of parameters to parametrize the remaining 
lattice artifacts. Figure [8] shows the fit results for some Z RGl utilizing different parameter 
sets {cfc}. We use the subtraction type (s) with g+ = g B . In that case the results from the 
complete one-loop subtraction (1) serve as reference values. 

Generally, we recognize that the resulting RGI renormalization factors do not vary 
significantly. Most fit results for Z RGl are located in a 1 % deviation band around the 
corresponding complete subtraction results (1). In addition, parametrizations (2) and (3) 
give almost identical fit results. This reflects, of course, the fact that our momenta are 
very close to the diagonal in the Brillouin zone. These restricted momentum sets might be 
the reason that even 'incomplete' hypercubic invariant sets (4, 5) can be used to obtain 
reasonable fits. For the final results we use the fit with all Cj^O which would be natural 
in the case of more off-diagonal momenta. 

In Figs. HI [10] and [11] we show the results for all operators using the parameter sets with 
all Ci compared to the results obtained by the subtraction scheme based on (1141) . 
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Figure 8: Z RGl for selected operators at /3 = 5.40 as a function of the parameters included 
in the fit ansatz ( l22i) . The used parameter combinations (partype) are: 
1: complete one-loop subtraction of lattice artifacts (114"]) 2: all q, 3: (ci,c 4 ,c 6 ) - 0(4) 
invariant, 4: (ci, c 2 , c 3 , c 4 , c 5 ) - (a 2 , a 4 )- hypercubic invariants, 5: (c 4 , c 5 , c 6 , C7, c 8 ) - (a 4 , a 6 )- 
hypercubic invariants. The horizontal borders of the shaded area show a 1% deviation from 
case 1. 
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6 Results for local and one-link operators and conclu- 
sions 



As a result of the preceding discussions we use subtraction type (s) (eq. (120]) ) with boosted 
coupling gs and the fitting formula ( l22i) with all q and h\ coefficients to determine the 
Z l . The final renormalization factors are collected in Table |2] using the two different 
r Aj^g values 0.700 and 0.789. This shows the influence of the choice of r A^jg (depending 
on the anomalous dimension of the operator). For the investigated operators and (3 values 
we found for the relative differences of the Z RGI 



5Z 



RGI 



7 RGI 7RGI 
r o a ms= - 700 r o A H g=0.789 



vRGI 

r A MS = 



: Q.700 



< 0.04 



(24) 



For comparison we collect in Table |3]the values for Z RG1 computed by means of fits with 
the ansatz ffl6l) to data where a complete one-loop subtraction of lattice artifacts (according 
to ffT4"j) with g+ = gs) has been performed. Note that here the errors are determined from 
the variation of the subtracted data between the scales /1 2 = 10, 20, 30 GeV 2 |4J. The 
reported renormalization factors are calculated for the values r /a given at the end of 
Section H] and, therefore, differ from those given in [31. The Z factors of the local operators 
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Op. 


r A-ms 


7RGI1 

£ 1/3=5.20 


7RGI1 

£ 1/3=5.25 


7RGI1 

£ 1/3=5.29 


7RGI1 

£ 1/3=5.40 


qS 


0.700 
0.789 


0.4530(34) 
0.4717(44) 


0.4475(33) 
0.4661(65) 


0.4451(32) 
0.4632(54) 


0.4414(30) 
0.4585(27) 


qV 


0.700 
0.789 


0.7163(26) 
0.7238(72) 


0.7253(26) 
0.7319(94) 


0.7308(25) 
0.7365(99) 


0.7451(24) 
0.7519(50) 


qA 


0.700 
0.789 


0.7460(41) 
0.7585(46) 


0.7543(40) 
0.7634(77) 


0.7590(39) 
0.7666(81) 


0.7731(37) 
0.7805(30) 


qT 


0.700 
0.789 


0.8906(43) 
0.8946(85) 


0.9036(42) 
0.9041(111) 


0.9108(41) 
0.9075(120) 


0.9319(39) 
0.9316(49) 


QV2,a 


0.700 
0.789 


1.4914(55) 
1.4635(108) 


1.5131(55) 
1.4776(112) 


1.5266(54) 
1.4926(90) 


1.5660(53) 
1.5397(58) 


QV2,b 


0.700 
0.789 


1.5061(37) 
1.4601(151) 


1.5218(37) 
1.4727(206) 


1.5329(36) 
1.4863(165) 


1.5534(35) 
1.5115(140) 



Table 2: Z RG1 values using the subtraction (s) with gs- 

in both tables agree within 1 %. The Z factors of the one-link operators differ at most by 
2%. 

Let us compare our results in Table [3] for the local vector current with Zy Gl obtained 
from an analysis of the proton electromagnetic form factor [19] following [20], which are 
listed in Table [H The numbers agree within less than 1 % with the numbers in Table [3] 
(r Aj^g = 0.700), supporting the complete one-loop subtraction as our reference point. 

From the present investigation we conclude: The alternatively proposed 'reduced' sub- 
traction algorithm can be used for the determination of the renormalization factors if the 
complete subtraction method is not available. Possible applications could be Z factors for 
Nf = 2 + 1 calculations with more complicated fermionic and gauge actions where one-loop 
results to order a 2 are available (for the fermionic SLiNC action with improved Symanzik 
gauge action see Ref. [TO]). 

In this study we have analyzed data sets with momenta close to the diagonal of the 
Brillouin zone. The one-loop a 2 contributions to the Z factors are completely general and 
can be used for arbitrary (also non-diagonal) momentum sets. Our ansatz ( |22l) allows to 
take into account the remaining artifacts after subtracting these one-loop a 2 terms. To get 
reasonable fit results the ratio (number of data points) / (number of fit parameters) has to 
be sufficiently large. 

As we pointed out the subtraction type is not unique. With (s) and (m) we tested two 
different types. The resulting fits do not give a clear preference for one of these. Even 
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Op. 


r A-ms 


7RGI1 

£ 1/3=5.20 


Z RGl /3=5.25 


7RGI1 

£ 1/3=5.29 


7RGI1 

£ 1/3=5.40 


qS 


0.700 
0.789 


0.4508(20) 
0.4620(85) 


0.44952(32) 
0.4603(60) 


0.44788(70) 
0.4585(61) 


0.4460(20) 
0.4560(48) 


qV 


0.700 
0.789 


0.7225(44) 
0.7219(53) 


0.7321(31) 
0.7316(41) 


0.7370(46) 
0.7364(55) 


0.7511(41) 
0.7506(50) 


qA 


0.700 
0.789 


0.7529(17) 
0.7530(14) 


0.76046(70) 
0.76054(48) 


0.76463(33) 
0.7647(14) 


0.77731(20) 
0.7774(10) 


qT 


0.700 
0.789 


0.9020(12) 
0.8948(40) 


0.91427(24) 
0.9072(32) 


0.9206(14) 
0.9137(48) 


0.94009(69) 
0.9333(38) 




0.700 
0.789 


1.5018(48) 
1.473(18) 


1.5190(64) 
1.490(14) 


1.5321(52) 
1.504(12) 


1.5681(29) 
1.540(14) 




0.700 
0.789 


1.5083(51) 
1.480(15) 


1.524(14) 
1.497(28) 


1.5362(92) 
1.509(23) 


1.5706(61) 
1.5436(69) 



Table 3: Z using a complete one-loop subtraction of lattice artifacts. 



7RGII 

£ 1/3=5.20 


7RGI1 

£ 1/3=5.25 


7RGI1 

6 1/3=5.29 


7RGI1 

6 1/3=5.40 


0.7296(4) 


0.7355(3) 


0.7401(2) 


0.7521(3) 



Table 4: Z RGl values for operator V from the proton electromagnetic form factor analysis. 

the additional choice for the coupling (g+ — g or g* — gs) does not lead to significantly 
different results. Therefore, our final choice (s) (eq. (120"]) with = gs) was supported 
by 'external' arguments: the improved behavior of the boosted perturbative series and the 
results obtained by complete one- loop subtraction |4]. 

We have shown that already the one-loop a 2 subtraction improves the behavior of the 
Z factors significantly: In the small p 2 region the contributions of the remaining lattice 
artifacts are smaller than the corresponding one-loop a 2 terms. As mentioned above, the 
accuracy to determine the Z factors is already at the 1 % level for local operators and at 
the 2 % level for operators with one covariant derivative compared to the complete one-loop 
subtraction of lattice artifacts. Additional systematic uncertainties are due to the choice 
of the r Ajyfg and r /a. 
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Appendix 

In this Appendix we show that the definition ([6]) leads to renormalization factors which 
are invariant under the hypercubic group H{4). 

We consider a multiplet of local quark- ant iquark operators Oi(x) (i = 1,2, ... ,d) in 
position space which transform according to 

OiW^SyWOjiR^x) (A.l) 

when 

ip(x) -> D(R) ijj{R- l x) , ip(x) ->■ i){R~ l x) D(RY (A.2) 

for all A" = 384 elements R of H(4). Here D(R) denotes the (unitary) spinor representation 
of H{4) (or 0(4)): 

D{R)^D{R) = R» vlv . (A.3) 

We assume that the operators Oi(x) have been chosen such that the d x d- matrices S(R) 
form a unitary irreducible representation of H(4). 

Denoting the unrenormalized vertex function at external momentum p of the operator 
Oi by Ti(p) we have 

d 

Ti(p) = Sij(R) D{R) F j (R~ 1 p) D(RY (A.4) 
j'=i 

for all R G -^(4), and analogously for the corresponding Born term Tf orn (p). Consequently 
we get 

d d 

J> [r,(p)r 4 ( P )t] = ]Ttr [r,(ifc)r,(i*p)t] . (a.5) 

i=l i=l 

Using the orthogonality relations for the matrix elements of irreducible representations one 
finds in addition 

1 d 

J> [r^r^t] = 2 6 ijJ2J2 tr [rk(R P )r k (R P y] , (A.6) 

R k=l R 

where the sum extends over all R G H(4). The same relations hold when one of the vertex 
functions or both are replaced by the corresponding Born terms, e.g., 

d d 

$> [r,(p)rf° r >)t] = [URp^iRpy] . (A.7) 

i=i i=i 

Therefore the renormalization condition 

z _i Etitrh^rfo^t] 
EtiM r f m b) r f ri W] 
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or, equivalently, 

z -i zs _ d E fi tr[r t (^)rf-(^)t] 

respects the hypercubic symmetry, i.e., writing more precisely Z = Z(p) we have Z(Rp) = 
Z{p) for all R 6 H(4), and all lattice artefacts in Z must be invariant under the hypercubic 
group. Of course, here it has been assumed that Z q {Rp) = Z g (p), as is the case for our 
definition Q of Z q . 
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